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We devise an ab initio formalism for the quantum dynamics of Auger decay by laser-dressed 
atoms which are inner-shell ionized by extreme ultraviolet (xuv) light. The optical dressing laser 
is assumed to be sufficiently weak such that ground-state electrons are neither excited nor ionized 
by it. However, the laser has a strong effect on continuum electrons which we describe in strong- 
field approximation with Volkov waves. The xuv light pulse has a low peak intensity and its 
interaction is treated as a one-photon process. The quantum dynamics of the inner-shell hole 
creation with subsequent Auger decay is given by equations of motion (EOMs). For this paper, the 
EOMs are simplified in terms of an essential-states model which is solved analytically and averaged 
over magnetic subshells. We apply our theory to the M4^5NiN2,3 Auger decay of a 3d hole in 
a krypton atom. The orbitals are approximated by scaled hydrogenic wave functions. A single 
attosecond pulse produces 3d vacancies which Auger decay in the presence of an 800 nm laser with 
an intensity of 10^^ W cm~^ . We compute the Auger electron spectrum and assess the convergence 
of the various quantities involved. 

PACS numbers; 32.80.Hd, 32.80.Fb, 32.80.Aa, 32.80.Qk 



I. INTRODUCTION 



The inner-shell ionization of atoms leads to a fasci- 
nating array of many-electron effects. Such vacancies 
decay on an ultrafast time scale by fluorescence or elec- 
tronic decay. Electronic decay refers to Auger decay [J 
and its special variant, Coster-Kronig decay 0- For 
low energy transitions in, e.g., light elements or high- 
lying inner shells in heavier elements, Auger decay is the 
dominant relaxation process For deep inner-shell 

holes, the ion relaxes in terms of cascades of fluores- 
cence and electronic decay processes. All atomic elec- 
trons are involved either directly in the photoionization 
with subsequent electronic or fluorescence decay or indi- 
rectly due to a rearrangement of the atomic electrons in 
the presence of newly formed holes, so-called core relax- 
ation 5]. Furthermore, at photon energies near the ion- 
ization threshold, the outgoing photoelectrons are slow 
and may interact significantly with the ionic remnant. 
If the subsequent electronic decay process is fast, out- 
going photo- and Auger electrons even repel each other 
appreciably. Thisphenomenon is dubbed post-collision 
interaction 1,1 Si. Generally, we will refer to elec- 
tronic decay and Auger decay and will not explicitly 
distinguish Coster-Kronig decay. They are fundamental 
processes that are pure manifestations of electron corre- 
lations. In this way, electronic decay processes are ideal 
for an investigation with the methods of attosecond sci- 
ence Jl, jLO, Uj, [13, [13 which aim to measure and control 
the motion of electrons on their natural time scale, which 
is the attosecond. 

Attosecond light was first used to measure the 
Auger decay time of 3d vacancies (Af shell) in 
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krypton atoms [ij, [l5|. Such vacancies undergo 
M4_57ViA^2,3 Auger decay which has been studied exper- 
imentally in the frequency domain [l^ [13, [3, [l^ ■ The 
time-domain study of Auger decay represents a semi- 
nal experiment in several ways. On the one hand, it 
demonstrates the power of the newly created attosecond 
methodology by comparing its results with existing data. 
On the other hand, the study of transient electron mo- 
tion with attosecond science has so far been restricted to 
mostly one-electron processes, e.g., Refs. [11,113, [HI and 
references therein. However, the most profound goal of 
attosecond science remains the study of electron correla- 
tions. 

Clearly, controlling a process on an attosecond time 
scale requires extreme ultraviolet (xuv) light for its short 
cycle period. Naturally, xuv light targets inner shells 
for which the photoabsorption cross section is highest 
at these wavelengths. The degree of control over inner- 
shell electron motion is potentially limited compared 
with valence electrons because of weak present-day at- 
tosecond light sources and the fact that postgeneration 
pulse shaping capabilities in the xuv and x-ray domains 
are severely limited, e.g., only an amplitude shaping can 
be accomplished so far with electromagnetically induced 
transparency for x rays j22|] . However, the attosecond 
light can be augmented by an additional optical laser, 
a so-called two-color problem. For moderate intensi- 
ties, the infiuence of the optical laser on electrons in the 
atomic ground state and hence the two-electron interac- 
tion among them can be neglected, e.g., noble gases can 
sustain very high electric fields before ionizing. For all 
elements, the impact of the optical laser on inner-shell 
electrons is negligible. For the optical laser to impact 
inner-shell electrons, its intensity needs to be so high that 
it would valence ionize the atom. Yet it can be used to 
structure the continuum, i.e., have an influence on liber- 
ated electrons, and in that way enable profound control 
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over electronic processes [IJ, [22| • A seminal experiment 
along the above-mentioned lines — but on a much slower 
(picosecond) time scale — is ultrafast laser control, using 
coherent excitation with short laser pulses, of the en- 
ergy and proximity of Rydberg electrons in an atom by 
Pisharody and Jones [23|. To demonstrate their ability 
to control electron dynamics, they showed that the au- 
toionization of the doubly-excited barium atoms is due 
to electron-electron collisions instead of a slow transfer 
of energy. Another notable work is the theoretical inves- 
tigation of an xuv pump-probe scheme for the study of 
the simultaneous two-electron emission in helium by Hu 
and Collins ^\. 

The time-domain measurement of the Auger decay 
of the 3d vacancy in krypton atoms has been ap- 
proached theoretically in two different ways [1^ [2^, [23] . 
First, Yakovlev and Scrinzi |2 5| devised a model to sup- 
port the experimental study |14l |. They made a simple 
rate equation ansatz to represent the 3d hole population. 
This population was used together with the Auger elec- 
tron wave function to replace the transition dipole ma- 
trix element in a formula to determine the laser-streaked 
photoelectron spectra of xuv-ionized atoms [2^. Second, 
Smirnova et al. [26| revisited the question of Auger decay 
by weakly laser-dressed atoms in terms of a fully coher- 
ent system of equations of motion (EOMs) formulation 
for an essential-states model [l^ using Hartree products. 
It is constructed to describe the Auger electron spectrum 
as simply as possible decoupling the EOMs in terms of 
a parametric decay width from Weisskopf-Wigner the- 
ory [33, HH [12] and solving the resulting system of equa- 
tions analytically. Based on the ideas in Refs. [25l|26l .l27| 
Zhao and Lin [33] and Wickenhauser et al. |34| stud- 
ied Fano resonances. Kazansky and Kabachnik [35l [36j 
developed an ab initio theory for the solution of the 
time-dependent Schrodinger equation for photoionization 
of inner atomic shells in terms of a Fano-Feshbach for- 
malism with short pulses that takes into account near- 
threshold effects. Finally, Smirnova et al. [13] applied 
the theory in Ref. [2^ to devise a scheme to use electron 
correlations to make attosecond measurements without 
attosecond pulses. 

Our study goes beyond previous work [H, [H, [13] and 
overcomes many of its restrictions. We develop a nonrel- 
ativistic multideterminantial ab initio formalism for the 
interaction of two-color light with atoms. We set out 
from the Hartree-Fock-Slater (HFS) approximation for 
the atomic orbitals. Such mean-field orbitals are typically 
a good representation to describe Auger decay [Q]. Us- 
ing general spin-singlet configuration-state functions, we 
derive EOMs. We treat the interaction with light semi- 
classically because, in contrast to a one-electron quan- 
tum electrodynamic formalism ^38,, .39i] , band width and 
pulse duration are treated more easily. Naturally, these 
play a decisive role in attosecond science. The general 
EOMs are subsequently simplified to an essential-states 
model [2^ and the equations are solved analytically in 
this special case considering, in contrast to the work in 



Ref. [26|, also the laser dressing of the photoelectrons. 
Furthermore, the laser dressing is, in our case, not re- 
quired to be weak, yet the intensity should remain below 
the excitation and ionization threshold of atomic ground- 
state electrons. We make a model for the atomic elec- 
tronic structure in terms of scaled hydrogenic wave func- 
tions [2^ [23]. Our formalism is a basis for the study 
of more complex situations in Auger decay and its con- 
trol. In forthcoming papers [i^, we will investigate the 
interference between Auger electrons from a twin xuv at- 
tosecond pulse. Further, one can examine what new av- 
enues for the control of Auger processes open up when 
one relaxes the assumption of an essential-states model 
to a multichannel treatment. 

The paper is structured as follows. In Sec. |TT1 we de- 
vise EOMs to treat the quantum dynamics of Auger de- 
cay on an ab initio level. Volkov waves are introduced 
in Sec. IIIII to describe the laser dressing. The EOMs 
are solved analytically for an essential-states model in 
Sec. [iVl We devote Sec. |V] to the determination of the 
dipole and two-electron matrix elements for our formal- 
ism. We then apply our theory to the laser-dressed Auger 
decay of krypton 3d vacancies; computational details are 
given in Sec. IVII and results are presented in Sec. IVIII 
Conclusions are drawn in Sec. IVIIII 

Our equations are formulated in atomic units [4l| . 
The Bohr radius 1 bohr = 1 ao is the unit of length and 
1 to represents the unit of time. The unit of energy 
is 1 hartree = lE/j. Generally, we use the indices h, i, 
j, m, . . . to denote occupied orbitals, a, b, c, d, ... for 
unoccupied orbitals, and p, q, r, s, . . . for orbitals which 
may be occupied or unoccupied. 



II. QUANTUM DYNAMICS OF 
PHOTOIONIZATION AND AUGER DECAY 



This section forms the core of our theory. In Sec. Ill Al 
we describe the schematic of XUV photoionization with 
subsequent Auger decay. The quantum mechanical foun- 
dation is laid out in terms of an ab initio description in 
Sees. Ill 51 III C[ and III D1 where we introduce the Hamil- 
tonian and the states involved. The full nonrelativistic 
formalism is simplified using an approximate Hamilto- 
nian in Sec. Ill El which comprises only those two-electron 
integrals which are essential for Auger decay. We correct 
for our omissions by adjusting the energies of the involved 
states appropriately. Finally, we use the time-dependent 
Schrodinger equation with the approximate Hamiltonian 
to formulate EOMs in Sec. lIIFI for the quantum dynamics 
of xuv absorption and Auger decay. 



A. Schematic of the processes 

Auger decay can be treated theoretically in various de- 
grees of sophistication \^ . We focus on a description 
in terms of a single isolated resonance (no decay cas- 
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FIG. 1: (Color online) Schematic of laser-dressed xuv pho- 
toionization of an inner-shell electron with a subsequent 
Auger decay. Only participating electrons are drawn. Photo- 
electron and Auger electron lines are influenced by a dressing 
laser and are observed in an electron spectrometer. Similar 
to Fig. la in Ref. 



cades). In the language of a full scattering process, Auger 
decay is a resonance in the double photoionization cross 
section . The scattering of an xuv photon 7 off an 

atom A leads to the formation of a dication A^^ and the 
emission of two electrons ep and e^, 

A + 7 — ^A+^+ep + e^. (1) 

Here, we ignore electron correlations among the ground- 
state electrons and in the cation and the dication. We 
also ignore the interaction between the outgoing elec- 
trons and the remaining ground-state electrons and the 
repulsion between the outgoing electrons, so-called post- 
collision interactions [B, @] . 

Assuming intermediate singly-ionized resonance states, 
we break up Eq. ^ into two separate processes: the 
photoionization of an atomic inner shell with subsequent 
Auger decay [Ijli^ which is depicted in Fig.[TJ To begin 
with, let us disregard laser dressing. A level scheme of 
the states participating in the photoionization and sub- 
sequent Auger decay are shown in Fig. [2] The atom 
is initially in the ground state with energy Eq- Then, 
it absorbs an xuv photon with an energy of wx- This 
leads to the formation of a singly inner-shell ionized 
cation with energy and the ejection of a photoelec- 
tron. With the Einstein relation, the peak of the energy 
distribution of the photoelectron spectrum — the nominal 
photoelectron energy flp — is found to be kp/2 = flp = 




► 

Interaction 

FIG. 2: Energy level scheme of the production of an inner- 
shell hole by xuv photon absorption and its subsequent Auger 
decay. The xuv photon energy is Wx, the atomic ground-state 
energy is -Eo, the singly-ionized state has an energy and 
the doubly- ionized state energy is E^^ . The photoelectron 
has a nominal kinetic energy of kp/2 = Qp and the Auger 
electron of k\/2 = 0,a. 



The inner-shell hole Auger decays; 
it is filled by a valence electron and the excess energy 
is transferred ultrafast by electron correlations to a sec- 
ond valence electron which is expelled. This gives rise to 
an Auger line in the electron spectrum at the nominal 
energy k\/2 = Qa = E^ - E^^ . Afterwards, the system 
is in a dicationic final state with energy E^^ . This ap- 
proximate mechanism of an ionization step with a fol- 
lowing electronic decay step is frequently referred to as 
two-step model of Auger decay d, |42| . Usually, this is a 
good approximation for atoms and molecules. However, 
in condensed matter, rearrangement processes take place 
which necessitate treating Auger decay as a one-step pro- 
cess [i,|42]. 

Next, we consider the influence of an optical dress- 
ing laser. Its intensity is assumed to be sufficiently 
low that it does not excite or ionize atomic ground- 
state electrons. However, the presence of the laser has 
a strong impact on the outgoing electrons and this has 
important consequences for the observed signal. We use 
Volkov waves [ij, IS] to describe continuum electrons 
and we exclude the possibility that the initial xuv ab- 
sorption does not ionize the atom but induces only an 
excitation to a Rydberg orbital. This means that ef- 
fects such as electromagnetic ally induced transparency 
for X rays [H, [H, HI, |45|, IH, |43| are not represented in 
our formalism. Always, the laser intensity is assumed 
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to be sufficiently low as not to modify the Auger decay 
rate noticeably with respect to the laser-free case. In 
other words, we neglect the impact of the laser on the 
energies i?o, E^, and E^'^ and the Auger decay-relevant 
two-electron matrix elements (three atomic ground-state 
orbitals and a continuum electron) because they are only 
weakly perturbed. 



B. Atomic electronic structure 



The contribution Hen in Eq. ([3]) allows for the fact 
that in the Hartree-Fock-Slater approximation ^ the 
electron-nuclear interaction and the electron-electron in- 
teraction are replaced by the one-electron central po- 
tential V^HFs(^) where the nucleus is at the origin of a 
spherical polar coordinate system and the electron is at 
radius r [38| . This replacement is reversed by the correc- 
tive term 



HCR = E ( I ^CH I <y5g ) [^pf ^gt + ' 



'pi ^i] • 



(5) 



We assume a Z-electron atom with a closed-shell spin- 
singlet ground state. The ground-state wave function 
is approximated by a Slater determinant of one-electron 
orbitals ^o{ficri, . . . ,fz crz)- Electron coordinates are 
given by fi and spin projection quantum numbers are 
given by for i e {1,. . . ,Z}. We use the formalism of 
second quantization where bp^. and b^^ are an annihilator 
and a creator of an electron in the spin orbital V'pct(^), 
respectively, which is the tensor product of a spatial or- 
bital ipj,(f) and a spinor of projection quantum num- 
ber a [41[. We define spin up cr = i = t and spin down 
cr = - i = J, and make the following ansatz for the atomic 
ground state, 



Z/2 

i'i'o)=n^it^iiio) . 



(2) 



using the vacuum state 1 ) [4l| . 

Let the full atomic electronic structure Hamiltonian 
be denoted by Hat- It consists of the kinetic energy 
of the electrons, the electron-nucleus attraction, and the 
electron-electron repulsion [4l| . We determine atomic or- 
bitals within the HFS approximation [4^|4^ to Hat- In 
a next step, we use the HFS orbitals to represent Hat- As 
the HFS approximation is typically good, Hat naturally 
decomposes into a part which is large and another which 
can be treated as a perturbation. We can rewrite if at 
as follows: 



Hat = Hups + Hqh + H^ 



(3) 



assuming a representation in terms of the atomic orbitals. 
We discuss these terms in the following paragraphs [4l| . 

The HFS Hamiltonian in Eq. ([3|) represented in terms 
of HFS orbitals reads 



H 



HFS 



>\b , 
pi pi 



(4) 



with 



hcH = -VHFs(f) 



Z 



(6) 



in addition to the explicit treatment of electron correla- 
tions in the next paragraph. There are no mixed terms 
involving a spin-up and a spin-down orbital in Eq. ([5|) 
because the interaction ([6]) does not depend on the spin. 
Electron correlations in Eq. ^ are described by 

it 

^ cr,cr',^,^'=lp.,p',q,q' 

_ _ (7) 

which is the only two-particle operator [32, |41| in the 
total Hamiltonian of the atom. The two-electron inte- 
grals in Eq. ([7]) are defined in terms of HFS spin or- 



bitals by Vpap'a'q^q'C = 

Coulomb interaction is 



Ifpqi^Pq'e)- The 



(8) 



where f and f' represent the coordinates of two electrons 
in the atom. 



C. Atom in xuv light 

We assume that the xuv light is linearly polarized 
along the direction ex and the wavelength is sufficiently 
large for the electric dipole approximation to be ade- 
quate jsol . Then, the Hamiltonian for the interaction 
of electrons with xuv light [50] reads 



Hx = ZDpq{t) [feJtV + ^PiM' 



(9) 



with the spin-independent orbital energies Sp. For nota- 
tional clarity, we assume a countably infinite set of final 
states in Eq. ([4|) and in this entire section, i.e., we assume 
a finite volume for box normalization of continuum wave 
functions [12] or a (finite-element) basis set expansion 
of the radi al p art of the atomic orbitals times spherical 
harmonics [3^. From Sec. IIV D] onward, we will use con- 
tinuum wave functions, namely, plane waves and Volkov 
waves fHlillil. 



using the matrix element Dpq(t) = ( (/5p ] /ix ] ) in terms 
of atomic orbitals of the spin-independent one-electron 
interaction with xuv light in length form, 

hx = r-Ex{t). (10) 

We use for the electric field of the xuv light 

Exit) = exit)exCOs{ujxt) , (11) 
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with the pulse envelope £x(i) and the angular fre- 
quency wx- With Eqs. PU|) and pT|) . we can rewrite 
the light-electron interaction as follows: 



have 



Dpg{t) = d 



pq 



(12) 



where the atomic dipolc matrix elements are given 
by dpg E ((/Pp|f•ex|v<^)• 
The complete Hamiltonian consists of the atomic 
electronic structure ([3]) and the interaction with the 
xuv light ([9]) and reads 



H = HxT + Hx 



(13) 



It will serve as the basis to treat the quantum dynamics 
of xuv absorption and subsequent Auger decay. 



D. Excited states 

Having formulated the Hamiltonian of the prob- 
lem ([T^ and having specified the ground state (H]), we 
need to incorporate singly- and doubly- excited states 
in our description to represent photoabsorption and 
Auger decay. We use a single configuration-state func- 
tion (CSF), which is a linear combination of singly- and 
doubly-excited determinants, to stand for a singly- and 
a doubly-excited state, respectively. Generally, singly- 
excited determinants are a satisfactory approximation to 
singly-excited wave functions which is the reason for the 
success of Koopmans' theorem |4]|. However, doubly- 
excited states are not so well represented by doubly- 
excited determinants because of hole-hole repulsion ef- 
fects and continuum electron interaction. To overcome 
this approximation, one needs to allow for configuration 
interaction 41]. 

Spin-singlet singly-excited states are represented 

by 



1 



x/2 



1 

7^ 



(14) 

Here, 6^^ creates a hole in the orbital h with spin projec- 
tion number a by destroying the electron and creates 
an electron in the orbital a with spin projection num- 
ber a. The hole orbital indices h which are taken into 
account form the set %; the indices for virtual (unoccu- 
pied) orbitals are a > -f (.^ is even because we consider 
only closed-shell atoms). After the second equals sign 
in Eg. (|14l). we introduce a concise determinantial nota- 
tion [41[. The bar over spatial orbital indices indicates 
a spin orbital with spin down; no bar refers to a spin 
orbital with spin up. 

There are five classes of doubly-excited spin-singlet 
configuration-state functions (Table 2.7 in Ref. [4l|). We 
focus on the two classes in which all four spatial orbitals 
are distinct. Using the concise notation of Eq. (|14p . we 



1 



[2|$^) + 2|$f ) 



'12 

+ l<ff^) + l<i>f ) 
1 



(15a) 



i*';) + i*f;) + i*-i') + i*!f )] -(isb) 



The pairs of orbital indices (i,j) of double vacancies 
with i < j, which are considered, constitute the set J^. 
With the restrictions a < b and a,b > ^ for the virtual 
orbital indices, we enumerate all distinct doubly-excited 
configurations. 

The energy of the ground state ^ is found with the 
electronic Hamiltonian ^ as follows: 

Zjl 

£;o = ( $0 I i^AT I $0 ) = 2 ^ + ( $0 I i^CH + i?cc I $0 ) • 

(16) 

The first term on the right-hand side is twice the 
sum of occupied Hartree-Fock-Slater orbital energies 
dubbed ^l|. The singly-excited states (|14l) have an 
energy of 

El = Hat I ) = £o-eh+ea+{ \ Hcr+H,c \ ) , 

(17) 

and the doubly-excited states have an energy of 



e: 



e: 



E. 



ah 



■ A -^ - n - V I I AT I \ ^''h 

I], A ij.B \ A^ij l-'-'ATlA^ij / 

= £o-e,-ej+ea + eb + {i^^^\HcH + Hee\i^--) ■ 

(18) 

We have E^^^^ = Efl^ because both doubly-excited CSFs 
consist of determinants with excitations from the same 
spatial orbitals into the same spatial orbitals, only the 
spinors change. As our nonrelativistic Hamiltonian ([3]) 
does not depend on spin, especially, it does not contain 
spin-orbit coupling, the energies are the same for all the 
excited determinants and thus also for the configuration- 
state functions A and B. 



E. Simplified Hamiltonian 

We have formulated a full ah initio description of the 
problem with a truncated excitation manifold [4l|. Our 
framework represents an ideal starting point for further 
simplifications. Eventually, it will be reduced to an 
essential-states model in Sec. IIVI that contains only the 
absolutely necessary energies and matrix elements to still 
describe the physics of the processes. We form a matrix 
representation of the Hamiltonian H [Eq. (|13p] in terms 
of the orthonormal basis, 



a,5>Z/2,a<fe,76{A,B}}. 



(19) 



We decompose the representation of H into an exactly 
solvable part Hq, given by the diagonal elements of 
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and a perturbation Hi , given by the ofT-diagonal elements 
of H. This is a so-called Epstein-Nesbet partitioning [5ll, 
Isl . [53I I. The exactly solvable diagonal part is written 
compactly in first quantization as 

Ho = i$o)i?o(<foi+ E E \'n)KCn\ 

hena>Z/2 

+ E E E \\^t-)E^,,{''^t!\, 

(i.j)ey^ a.b>z/2 ^e{A,B} 

a<b 

(20) 

where we use the energies from Eqs. ((T6)) . (fT7| . and p8|) . 

These energies are well suited to be treated as (exper- 
imental) parameters [see Sec. IVIj . Here, H^ makes no 
contribution to the diagonal matrix elements of H be- 
cause it consists only of off-diagonal elements. 

The perturbation is given by Hi. We do not use all 
off-diagonal elements of H for Hi and, additionally, we 
make approximations to the ones we retain. We use 

Hi = z z [\'i>o){'i>o\Hx\'n){'n\ 

heH a>Z/2 

+ \'n)Cn\Hx\^o)(<^o\] 

h^H o.,b>ZI1 -yefA.B} 

(21) 

Here, Hufs [Eq. Ill] makes no contribution because it 
consists only of diagonal elements. The impact of Hen 
[Eq. (O] is neglected totally and we consider the energies 
in Eq. (jSO]) to be parameters which shall compensate for 
this and the other omissions that influence the energies of 
the involved states. Similarly, electron correlations Hcc 
are only included when they cause transitions between 
singly- and doubly-excited states, i.e., they are taken into 
account when they are responsible for Auger decay which 
cannot be understood in a mean-field picture. In princi- 
ple, the neglected matrix elements can be incorporated 
allowing one to carry out a full ab initio treatment of the 
problem. 

F. Equation of motion formulation of 
photoionization and Auger decay 

To describe a photoionization process with subsequent 
Auger decay, we solve the time-dependent Schrodinger 
equation 

i?|vE',0 = i^l*,t) (22) 

for an atom in xuv light. In terms of the states in the 
basis B [Eq. (fT9|) ]. a general state ket (or wave packet) is 
given by 

= co(Oe-'^°*|*o)+ E c^(i)e-'^^*r*J'J 

+ E E c£SV'^-^*i^$^), 

a,b>Z/2 
-ye{A,B} a<b 

(23) 



which we insert into the time-dependent Schrodinger 
equation ([^ . Exploiting Ho\(l)) = E^\<j)) ior states | </> ) e 
B with energies £'</,, we arrive at the EOMs for the expan- 
sion coefficients co{t), c^(t), and cfj j{t) by projecting 
on {(pl for all 1 ) 6 i3. The atom is initially in the ground 
state which implies the initial conditions co(0) = 1 and 
<(0)=c°,^^(0) = 0. 

We get the first EOM for ( 1 = ( $0 I which represents 
the rate of change of the ground-state amplitude, 

co(i) = -ix/2E E A.a(t)c'(^"-^^"^*c^(t). (24) 

hen a>Z/2 

We consider here the weak absorption limit, i.e., co(t) ~ 
1 for all times. The rate of change will, nevertheless, 
prove highly beneficial in determining the cross section 
in Sec. IIV El and the photoelectron spectrum in Sec. lIVFl 
The second EOM results from |0) = 1^$^) and de- 
scribes the inner-shell hole amplitude for h eH, 




(25) 

The first term on the right-hand side of the equation rep- 
resents hole production due to absorption of xuv light; 
the second term describes the loss of hole amplitude 
caused by Auger decay. In this equation, the two-electron 
matrix element in terms of spatial (e.g., Hartree-Fock- 
Slater) atomic orbitals is denoted by 

Vpp'qq' = ( ifp ipp' I hcc I Vq Vq' ) ■ (26) 

Further, we define the antisymmetrized two-electron ma- 
trix element Vppi[qqi^ = Vpp'qq' - Vpp'q'q and the sym- 
metrized two-electron matrix element Vppi^qqfj = Vpp'qqi + 
Vpp'q'q which consist of a direct matrix element Vpp'qq' and 
an exchange matrix element Vpp'q'q. 

The third EOM is obtained setting ((/)| = {\<^°;-\; it 
describes the Auger decay amplitude 

= i\/f E^'.''fo]e'<''^'*-''^°'*4(0,(27a) 

i^l^it) = i^ E«'«;>to}e'^''^'^-'^^^*c^(t).(27b) 

We can reduce the EOMs (|^ and (P?)) further by ignor- 
ing the electron exchange matrix element v^bji and ex- 
ploiting the fact that the energies of the doubly-excited 
states (UHl) are the same, Ef^ = Ef^j^ = Ef^^. Then, 
the second EOM (j25p for the inner-shell hole amplitude 
simplifies to 

+ i2V^ E E e'(^^0*c^(t), 

(28) 
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for h e H with the definition c°''(i) = c°';B(t). Again 
ignoring electron exchange and using Efj', leads us to 
the relation cfh{t) = V5cf'j'j^{t) between Eqs. (I?7a|) and 



([m]) : thereby, we use that both c'^l^it) and b(0 

van- 
ish initially. Then, we need to retain only the simplified 
Eq. (j27bp of the two EOMs for the Auger decay ampli- 
tude (1271) yielding 



•ab 



he-H 



(29) 



Equations (|28p and (l29l) constitute a linear system of dif- 
ferential equations which contains all phase information 
and thus describes interference effects. 



where the amplitude is Al,o = --^lgl- The laser phase 
at t = can be specified using At. Then, the laser electric 
field follows from 

^L(t) = -^^%^ = ^L,ocosK(i + At)), (32) 
at 

where the electric-field amplitude is -El.o = '^l -^l gl [see 
Eq. (pTl) ]. We assume free fields, i.e^a vanishing scalar 
potential, and the Coulomb gauge |4J, |50] ■ 

The Hamiltonian of a free electron in a laser field in 
velocity form [3, [s^l is 



(33) 



III. LASER DRESSING 

In Secini we devised a formalism to describe the quan- 
tum dynamics of the photoionization of the inner shell of 
an atom by xuv light and the subsequent Auger decay. 
Here we expand our formalism to include an additional 
optical dressing laser of moderate intensity. The impact 
of the laser on ground-state electrons is neglected and 
only the modification of the continuum wave functions of 
the photo- and the Auger electron is considered. 

To begin with, we simplify the manifold of virtual 
states by replacing it by free-electron wave functions, the 
momentum normalized plane waves [s^l, 



(27r)3/2 



(30) 



This substitution expficitly excludes Rydberg states. It 
is justified by the fact that we are only concerned with 
continuum electrons of sizable kinetic energy. A conse- 
quence of our replacement is that the new continuum 
wave functions are no longer strictly orthogonal to the 
bound-state wave functions because both sorts^ of wave 
functions stem from different Hamiltonians: ^ = 

and if HPS [Eq- O] 113 • Note that we exploited strict 
orthogonality in the derivations of Sec. HH 

When we additionally consider a laser field, our re- 
placement of continuum wave functions becomes known 
as strong-field approximation [44] in which the influ- 
ence of the Coulomb potential on continuum states is 
neglected. In other words, laser dressing can be incorpo- 
rated easily into our treatment by replacing the laser-free 
continuum functions (pp ]:{r,t) by Volkov waves [43l |4^ 
[see Eq. dSH) below]. In doing so, electrons in the atomic 
ground-state orbitals, however, are considered to be un- 
influenced by the laser. 

The laser pulse is assumed to be long with respect to 
all other time scales in this paper and is taken to be 
monochromatic and continuous wave. Let the laser radi- 
ation of angular frequency o^l be linearly polarized with 
the polarization vector cl. The vector potential is 



^L(i) = Al,o sin(wL(i + At)) , 



(31) 



using the electron momentum operator p = -iV. The 
solution of the time-dependent Schrodinger equation 
with Hamiltonian (1551) reads 



1 



^v,^-(nO=(^en—v ^^^^^^ (34) 

and is called a Volkov wave [i^ H, [Hi]. The Volkov 
phase is given by 



t 

,(fc,t) = iy" [k + Ai^{t')fe''''dt'. (35) 



The factor e''* for 77 > ensures the convergence of 
the integral. After the integral has been performed, the 
limit 77 ^ 0^ is taken. The Volkov phase vanishes at -00. 
Inserting Eq. dST]) into Eq. dSHD, we find 



(fc,<) = (^ + C/pjt-fc-aL cos(tJL(i + At)) 



Up 

2WL 



sin(2wL(i + At)) 



(36) 



where the ponderomotive potential is 
Up = — — = — alh,o , 



(37) 



with the electric-field amphtude |i?L,ol = v^StTo"/]^ for a 
laser with intensity /l,o- The fine-structure constant is a. 
During a laser cycle, the maximum classical excursion 
from the origin of a free electron is given by [4^ ] 



5 . 



(38) 



The exponential of the Volkov phase e-i*v(fc,t)^ 
can be expanded using the generating function of the 
generalized Bessel functions J,n{u, v) (44. Is^ which reads 



^-i {u cos 0+v sin(2 6)) 



E i-ir e'"'' Jm{u,v) , (39) 
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by setting = 6* - f in Eq. (10) in Ref. The Jm{u,v) 
can be evaluated in terms of the ordinary Bessel func- 
tions Jm-2n{u) and Jn{v) usiug [il. [55II 



Jm{u,v)= ^ J7n-2n{u) Jn{v) 



(40) 



Another connection to ordinary Bessel functions repre- 
sents the relation Jm{u) = Jm{u,0) with which Eg. ((59)) 
reduces to the familiar Jacobi- Anger expansion [56j. Fi- 
nally, the temporal phase factor reduces to 



E 



^-im (wl (t+At)-f ) 



(41) 



IV. ESSENTIAL-STATES MODEL FOR 
LASER-DRESSED PHOTOIONIZATION AND 
AUGER DECAY 

This section is devoted to a solution of the EOMs from 
Sec. Ill FI for an essential-states model [HjIl^liSl which is 
a fairly drastic approximation to the model Hamiltonian 
of Sec. Ill El However, it retains the essential physics of 
the problem. Namely, we include only the states from 
three magnetic subshells of the occupied orbital mani- 
fold and form essential states from them by averaging 
the dipole and two-electron matrix elements in Sec. IIV Al 
over the subshells. Subsequently, we decouple the sys- 
tem of EOMs in Sec. IIVBI using perturbation theory. 
The resulting equations are adapted to account for laser 
dressing following Sec. IHII and are solved analytically. 
We obtain the hole-state amplitude in Sec. IIV Cl which is 
used to find the Auger electron spectrum in Sec. IIVDI 
Finally, we derive the xuv absorption cross section in 
Sec. IIV El and the photoelectron spectrum in Sec. lIVFl 



A. Matrix elements and energies 

To begin with, let us simplify the problem and the no- 
tation. We represent photo- and Auger electrons by plane 
waves (|30p with momentum vectors fcp and ^a, respec- 
tively (laser dressing is not treated at this point). More- 
over, we make the following replacements: a fcp and 
& ^ fcA in our previous equations. Summations over a 
and b become integrals over fcp and fcA- As /i, i, and 
j refer to individual orbitals, their use is not meaning- 
ful anymore in our model context and are eliminated or 
adequately substituted as detailed in the following. 

To construct the matrix elements of the essential-states 
model, we note that the set of hole orbital indices T-L 
refers to orbitals from a single magnetic subshell. The 
pairs (j, j) e refer to orbital i from one magnetic sub- 
shell and orbital j from another magnetic subshell. The 
number of single hole states is #7^ and the number of 



double hole states is The xuv interaction matrix 

element is taken to be 



d{kp) = Qd 



1 



while the Auger decay matrix element is 



w(fcA) 



A 



— - — E E 

ifH #JF 



h fcA ijl 



(42) 



(43) 



Here, > is the strength of the dipole matrix ele- 
ment and Qy > is the strength of the two-electron ma- 
trix element. Both strengths will be determined later in 
Eqs. ([701) and respectively, based on (experimental) 
parameters. We chose to use the rms value to form av- 
erage matrix elements because in the following equations 
frequently the modulus squared of the matrix elements 
is used. 

Within the scope of our essential-states model, the en- 
ergies of the states of Sec. Ill Dl and Fig. [5] are as fol- 
lows. The ground-state energy ^6]) is Eq = Sq, neglecting 
the influence of Hqh + H^c- The energy of singly-excited 
states (|17l) is decomposed into the kinetic energy of the 

photoelectron and the energy of the cation E . It be- 

^2 

comes Ef^^ ^ + E^ with E^ = So-st (the orbital indices 
in V. denote orbitals from the same subshell and thus Sh is 
the same for all h e T-L). For doubly-excited states (|18p . 



we set E^^ ^ + E^+ with = So - e, - ej being 

the energy of the dication [again Si and ej are the same 
for all (i,j) e !F]. With these definitions, we find for the 
single ionization potential = E^ - Eq = -Sh and for the 
double ionization potential I^^ = E^^ - Eq = -Si - Sj [see 
also Fig. [5]. The nominal photoelectron energy from the 
stationary-state energy level scheme is ilp = Eq- E'^ + ujx- 
Likewise, JIa = E'^ - E^^ is the nominal Auger electron 
energy. 

Under these assumptions, we find for the inner-shell 
hole amplitude (pS)) — in which we replaced the hole in- 
dex h, by the subscript P and the double index ij by the 
subscript A — the expression 



(^-Op)* 



— d(fcp)ex(t)e'^- 



-i{-^-QA)t -fcp fcA 



(t) d^fcA 



= w{kp,t) 
The Auger decay amplitude ((29|) is 



(44) 



gfcpfcA(^) = ^«(fcA)e'(^-"*)*4^(0 . (45) 
V 2 

The EOMs of the essential-states model are formed by 
Eqs. (gH) and (gS]). They are very similar to Eqs. (10) 
and (11) in Ref. ^,|52|. 
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B. Perturbative decoupling 



Despite our considerable simplifications in Sec. IIV Al 
Eqs. (|44|) and (|45l) still form a linear system of cou- 
pled integrodifferential equations. The coupling stems 
from the second term w{kp,t) on the right-hand side 
of Eq. (UH). It describes the Auger decay of inner-shell 
holes and can be approximated in terms of second-order 
time-dependent perturbation theory (Weisskopf-Wigner 
theory) [s^, [Ml, [s^l . This treatment allows us to decouple 
the differential equation (|44p by eliminating the depen- 
dence on c'^ '^^ (t) in terms of a decay width F and an 
energy shift Ar as follows: 



w{kp,t) = [-iAn--)c^{t) . (46) 
The energy shift of the resonance state follows from 



(jj)6.?-7e{A,B} ^3 



l^kp V |2 



j-ikp T-ifcp kA 



■d-'/cA 



(47) 

where Pr indicates that the principle value of the integral 
has to be taken. The result was obtained by neglecting 

electron exchange and using K^f = E^fl^ = E^f^^. 
With the same assumptions, the decay width becomes 



(ij)€jc-^6{A,B} 



-2^ E /2K,,.,P'5(£;^-4-^--)d3fcA. 

m 

As we assume the two-step model of Auger decay [a, 143] , 
Th is independent of the photoelectron momentum fcp. 
In Eqs. (|47p and (|48p . we mark explicitly the dependence 
on the hole orbital h e Ti. As all /i e "H are from a 
single magnetic subshell, the Ap ji and the Th agree for 
all h e Ji. Therefore, we may drop the dependence on h 
in what follows. 

In the derivation of Eqs. (gH]), (gT]), and (gS]), we im- 
plicitly assume that c'^ (t) varies only slightly on time 
intervals [t - r; t] for t e [-oo; oo] and a small r > with 
respect to all time scales in the problem [g^, • Other- 
wise Eq. (|46p would not be meaningful. We can integrate 
Eq. (gS]) formally, 



L L— I 



t-T 
1 

7^ 



i(-^-r2A-ii;)t 



(49) 



g2 

-^-VLx-\ri 



Here, e''* with > ensures the initial condi- 
tion c^^'^*(-oo) = and the convergence of the integral 
where ^ 0^ is performed after the integration. The 



last equahty follows from / t'}^''^{t') dt' 



=fep kA 



T;kr kA 



{t -t), where the term for t - t cancels the second 
integral. 

Our result for the Auger amplitude (|49p is inserted 
into the expression for 'w{kp,t) [the second term on the 
right-hand side of Eq. (gl])] yielding 



.(fcp,t)=2i f iyA)pe''* ct(t) d^fcA 

-I k. ^ . 



^ - f^A - i?? 



With the decomposition [5£ 



— ^— = Pr - + m6{x) 

X - IT] X 



(50) 



(51) 



we obtain Eqs. (gll), (gT]), and (g8]) after dropping the 
subscript h, eliminating the sum over final states and 
replacing the energies and two-electron matrix elements 
in Eqs. (gT]) and (gS]). In detail, we find 



and 



T = 27T J 2\vCkA)\'S{nA-^)d!'kA 



(52) 



(53) 



A suitable value for the strength of the Auger decay 
matrix element (|43p can be obtained from Eq. ([55]) via 



(54) 



provided the decay width Ppar is taken to be an (experi- 
mental) parameter. 



C. Hole-state amplitude with laser dressing 



In Sees. IIV Al and IIVBI we disregarded laser dress- 
ing and focused on the EOMs with xuv light only. 
In the framework of Sec. IIIII we can easily incorpo- 
rate laser dressing in the strong-field approximation [4j| 
into our equations; the only change in our EOMs (|44p 

and (|45p concerns the time-dependent phase factors i^ t 
which need to be replaced by Volkov phases i^y{kp,t) 
[Eq. ([M]) ]. Using relation (gSjl. we decouple the hole 
amplitude from the Auger decay amplitude and recast 
Eq. (gH) into 



c^P(i) = -Ld(fcp)£x(t)e'*^(^^'*)e-'"P* 
v2 



+ (-iAR--)4-(t), 



(55) 
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assuming that the Auger decay is uninfluenced by the 
laser and thus the second-order energy shift Ar [Eq. ([5^ ] 
and the Auger decay rate F [Eq. ([53]) ] are meaning- 
ful. The first-order ordinary differential equation ([5l 
is solved analytically ^] yielding 



— oo 

^ gi*v(fep,t) gi(AR-i|-np)t' ^ 



To solve the time integration in Eq. (j56p . we expand 
the Volkov phase as in Eq. (|4ip and insert the inverse 
Fourier transform of the envelope of the xuv pulse pT|) , 



oo 

ZTT J 



(57) 



We obtain for the laser-dressed hole amplitude the 
expression 



following closed-form expression: 



-1 

47r 



d(fcp)i;(fcA) E e'™("-^*-t) 



oo 



Up ' 



ex(w) 



+ TO cjl + Up - rip + Ar - w - i- 



gi (-f-+mii;L + C^p-i2p-i"2A-'^) *' 

— oo 

(59) 

We are only interested in the Auger electron spectrum 
after the xuv pulse is over and the induced hole ampli- 
tude has decayed. Therefore, after expanding the Volkov 
phase (|4ip . we let i oo and simplify the time integra- 
tion in Eq. ([55]) by observing that 



S{(jj - to 



oo 

') = — [ e' di' 
' 2ti J 



(60) 



is a representation for Dirac's 5 distribution [33,l56[. Fi- 
nally, we obtain the laser-dressed Auger amplitude, 



'{t) 



2\/27r 



J(fcp) E e 



i 771 (ljl Ai--J) 



-Alp k a 



Jm^aL ■ kp 

ixi^) e 



Up 



) 



-+m uJi^ + Up -uj) t 



+ mwL + Up - flp + Ar ■ 



- (58) 

For moderate laser intensities, we have r* 0. 

Then, the generalized Bessel functions go over into 
ordinary Bessel functions [H, [5^. Further, the 
limit lim J^^u^v) = S„i q exists which completely re- 

moves the dependence of the equation on the laser for 
vanishing intensity. With this approximation, our ex- 
pression (158]) goes over into Smirnova et al.'s [H, 
Eq. (15). 



D. Laser-dressed Auger electron spectrum 



The amplitude to observe an Auger electron with mo- 
mentum fcA for a photoelectron with momentum fcp at 
time t is found by integrating Eq. (|45p from -oo to t. Be- 
forehand, expression (|45|) needs to be adapted for laser 

dressing by replacing i-^t with i<i>v(fcAji) after which 
we insert the hole-state amplitude (|58p. We obtain the 



(oo) = -d{kp)v{kA) E e 



i (m+n) ^i-'^) 



^ 2 wl ' ^ 2 cjl ' 

,p p . 

5'[-^ + TOtJL + C/p,-^+r7.WL + C/p) , 



with the line shape function 



(61) 



or \ £x(^p+^a-^p-»a) 
S{ujp,UA) = T — , (62) 

UJA - "A - Ar + 1 2 



which depends only on the absolute values of the mo- 
menta fcp and fcA- Formula (|6ip goes over into Eq. (18) 
in Ref. [2^ — apart from a factor ^ in our expression — by 
setting J7p = 0, replacing the generalized Bessel functions 
by ordinary ones and removing the dependence on the 
laser for the photoelectrons by using Jm(0) = Sm,o- 

We are now in a position to determine the laser-dressed 
Auger electron spectrum where we consider the case that 
the photoelectron is not observed. Therefore, we inte- 
grate the probability density jc^^ ''*(oo)p [Eqs. dS]) and 
((62|) ] over all possible photoelectron momenta to elimi- 
nate this degree of freedom. This yields for the prob- 
ability density [3^ to observe an Auger electron with 
momentum vector fcA, 



PA(fcA) = / |ci"^M°°)l'd^fcp 

R3 



(63) 
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E. xuv absorption cross section of laser-dressed 
atoms 

The probability of finding an atom in the ground state 
is given in terms of the ground-state amplitude co(t) in 
the wave packet (pS)) by 



Vo{t) = \co{t)\'' = c;{t)co{t) 



(64) 



Consequently, the negative of the xuv-absorption 
rate [3^ is 

« t*{t) + to{t) (65) 
= 2 Reco(i) = 2 lui[ito{t)] . 



The center line follows from the weakness of XUV absorp- 
tion, i.e., co(i) ~ 1 for all t. The rate of change of the 
ground-state amplitude follows from the first EOM ((24| ; 
adapted for the essential-states model with laser dressing, 
it reads 



(66) 

Inserting this EOM into Eq. ([55]) . expanding the Volkov 
phase factor using Eq. (|4ip . and inserting the laser- 
dressed hole-state amplitude ([55)) . we obtain the rate 



rx(<) = -2 Im [ico(t)] = — ex(i) / JmfaL-fcp,- ) Jnf (3l ■ fcp, ) 

27r m.n=-oo ^ 2lJl'^ ^ 2 LJL 



|J(fcp)P / 



Im 



£x(w) e 



ioj t i (ri-m) [ojl (t+At)--|-] 



^2 

^ + nWL + ?7p - f^p + Ar - W - i^ 



(67) 



dijj d fcp 



The absorption rate ([67|l in conjunction with the flux Jx = at the xuv (peak) intensity /x o with photon 
energy wx allows one to obtain the xuv photoabsorption cross section ^60] via 

Fx 



(68) 



Note that for a continuous- wave approximation of monochromatic radiation, we have for the xuv field strength ex(w) = 
ex.0 a-nd ex(i) = ^'X.o = ^87ra/x,o = const for all t. The relation between time- and frequency-domain field 
amplitudes follows from Eq. (j57p and is i?x,o = 57 £x,o- For large terms with m + n oscillate rapidly in the absorption 
rate (|67p . We discard these terms and retain only the constant terms with m = n. This leads to the expression 



(t(wx) = Siraujx 



00 ^ 
''=-°° R3 



(3l ■ fcp 



£/p 

2WL 



)^ Im[ 



g2 

i?+ + ^ + m cjl + [/p - i?o + Ar, - wx - i-ij 



(69) 



by expanding flp = Eq - + wx • For vanishing laser in- 
tensity, the structure of this equations becomes the same 
as from Eq. (40) in Ref. [11]. There, however, Ar was 
not accounted for. 



Similarly to Eq. ([54]) , a suitable value for the strength 
of the dipole matrix element (|^^ can be obtained from 
Eq. (I69p . omitting laser dressing, via 



Qd 



''par 



(l-^par) 
\ Cr(Wpar)|Qd = ] 



(70) 



F. Laser-dressed photoelectron spectrum 

We determined the rate with which xuv light is ab- 
sorbed TxCO [Eq. ([65])] in Sec. ITVEl The rate was de- 
rived under the premise of weak xuv absorption which 
allowed us to approximate the ground-state amplitude 
by co(t) ~ 1 for all t. Therefore, |co(i)P cannot be used 
to obtain the probability with which photoelectrons are 
ejected. Instead, we need to integrate the rate 

t t 

rp{t) = J rx(t') dt' = -2 J Im [i go(t')] dt' . (71) 



provided the cross section crpar('jJpar) at an energy Wpar 
in the range of energies of interest or close to the range 
is taken as a (experimental) parameter. 



In expression ([71]) . we insert the ground-state amplitude 
rate of change (j66|) and the laser-dressed hole-state am- 
plitude (|58p and expand the Volkov phase factor (|4ip to 
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obtain the probability 



t oo 

'Ppit) = 1^ / Jm(aL-kp,- )jJ&L-kp,- ) / ex{t') / Im 

271" „.^_oo ^ 2wl^ ^ 2uji^^ J J 



\d(kp 

X 4-^ : 



|2 -ii^t i(n-m)[wL(t+At)-f] 



+ n UJ]^ + Up - rip + Ap; - UJ - 



di' d^fcp 



ex(t^) 



(72) 



Letting i ^ oo, replacing the real- valued £x{t') by the complex conjugate of Eq. (|57|) . and omitting the integration 
over fcp, we find the probability density for photoelectron ejection by the xuv pulse. The time integration yields a 
6 distribution (|60p of the form S{lu' - uj + {n - m) cjl)- Replacing rh = n - m, we arrive at the probability density 



1 

2^ 



Pp(fcp) = -^|d(fcp)|' 



y. y, Jn-rh(aL-kp,^^)jn(aL-kp,-^^) 



Up 



7n=-oo n=— c 



2UJ-L' 



2WL' 



X Im 



2 „ 

f- + n ui^ + Up - ftp + Ap - u! - i=T 



(73) 



r 



The photoelectron spectrum depends on the Fourier 
transform of the xuv field envelope at lo and at w-mwL- 
This functional dependence indicates interference effects 
between channels with a different number of laser pho- 
tons provided that the xuv field envelope has sufficient 
width. 



V. ELECTRONIC STRUCTURE 

The theory of Sees. [TTl [Till and llVl treated the electronic 
structure of an atom as an abstract quantity which was 
represented by the orbital energies in Hhfs [Eq. (51)], 
the one-electron matrix elements in Hqh [Eq. ©J, the 
dipole matrix elements in Hx [Eq. and the two- 

electron matrix elements in Hee [Eq. Programs 
exist to carry out the Hartree-Fock-Slater approxima- 
tion m, |4§| and compute the required one- and two- 
electron matrix elements. To evaluate the essential-states 
model of Sec. llVi however, we use a much simpler model 
approach in terms of scaled hydrogenic functions for the 
atomic orbitals. This treatment follows Refs. psi. [26l. \2T\ . 
The parametrization of the model corrects to a large ex- 
tend for inaccuracies in the orbital energies and the ma- 
trix elements. If the results of the essential-states model 
depended sensitively on the electronic structure, then 
due to the substantial simplifications made, its physical 
predictions would be untrustworthy. Despite the use of 
approximate orbitals, the equations derived in this sec- 
tion are completely general and an ab initio evaluation 
in terms of Hartree-Fock-Slater orbitals is feasible. 

We use hydrogenic wave functions to model the 
spatial atomic orbitals in spherical polar coordi- 
nates Lpi{r,-d,Lp) = Rn^i^{r) Yi^jni{'&,v) lH- Here, Ui, 



li, and rrii are the principal, orbital angular momen- 
tum, and magnetic quantum number, respectively, of or- 
bital i 6 {1, . . . , Z/2}. The radial part is Rmhir) and the 
angular dependence is described by spherical harmon- 
ics Yi-miiiS,^)- We scale the hydrogenic wave functions 
such that their energy = matches the energy £i 

of the corresponding orbital in the chosen atom [2^, [23| . 
For this purpose we use an effective charge 

Zeff,^ = n,^/-2e^ . (74) 

This scaling also adjusts the spatial extend of the orbital 
appropriately. 



A. Dipole matrix elements 

The dipole matrix elements for xuv absorption (|12p 
are given by the promoted wave function in momentum 
space [3^ |6]|. 

'^SP. = 7^ / ^-''-'-Mr)d'r, (75) 

with f ■ ex = z. We use an atomic orbital <^/i(r) for 
the vacancy created by photoionization and the spatial 
part </Jpg(f, 0) of a plane wave (|30l) . This comprises 
also the case when laser dressing is considered because 
the laser dressing manifests itself exclusively [13] in the 
time-dependent Volkov phase ([551) . We use the Rayleigh 
expansion |62|] of the plane waves in Eq. ([751). 

00 I 

e''-' = 4^E E i'Yi:^{dk,^k)Mkr)Yi„.{^,^) . (76) 

1=0 m=-l 
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The directions of k and r are specified by the polar an- spherical Bessel function [5g|. We arrive at the dipole 
gles i^fc, tpk and 1}, ip, respectively. Here, ji denotes a matrix element (j75p in spherical polar coordinates, 

I 



(77) 



Corresponding to orbital /i, we have the principal n^, orbital angular momentum 1^, and magnetic quantum 
numbers. The angular integral is 



4tt 



(2ii + l)(2;2 + l) 
47r 



(2;3 + i) C{li,l2,h;mi, 7712,1713) C(/i,/2,^3;0,0,0) , 

I 



(78) 



where C(h , h\ wi , TO2, m^) is a Clebsch-Gordan coef- 
ficient [64l. The integral restricts the accessible angular 
momenta and magnetic quantum numbers in the pho- 
toionization process. The radial dipole matrix elements 
are 

00 

^Si.(^P) = / Mkpr)r'R„,i,{r) dr , (79) 


in terms of the radial part (r) of the atomic or- 

bital h. 



B. Auger transition matrix elements 

Auger decay is mediated by the two-electron matrix 
element, 



— 

''^h fcA ij ~ 



CkA) = j^fe-^'^'^^{f)d'r. (80) 



Here, x^^(r) and Xh{kA) are the configuration space 
and momentum space Auger electron wave functions, re- 
spectively [2^ [s^l . With the two-electron repulsion hcc 
[Eq. ([8])], the configuration space Auger electron wave 
functions reads 

00 

Xh{r) = <fj{r) J J (pl{r',Q.')hcc g}i{r' ,D.')r''^ dr' dVl' . 



4ir 



(81) 

To simplify h^c, we replace it by the Laplace expan- 
sion |63l |. 



1 00 A I 

S-7j = E577T«(M-')E/i;.K,.')«„{<>..) 



with the decomposition 



(82) 



Qi{ry) = —e{r-r') + — e{r'-r), (83) 



for the radial dependence, where 6 is the Heaviside step 
function with 0(0) = The wave function of the Auger 
electron is in spherical polar coordinates, 



47r 



L with A 



ihi'i) 2Z + 1 



X y{lh,lJi'i'm,h,'mi - mh^rrii) 



(84) 



The symbol A{lh lli) represents the triangular condition 
for which Clebsch-Gordan coefficients do not vanish [62| • 
The radial dependence in Eq. (fM]) is expressed by 



00 

<'Ln,iJr) = f R^,^,,ir')Mr,r')R,,^Ur')r"dr'. 


(85) 

We would like to calculate the momentum space 
representation (I80p of the Auger electron wave func- 
tion (gll) [13,161]. The plane wave is expanded in terms 
of spherical Bessel functions (|76l) . 



1 



4^2^ y 

^ 9 7+1 

i with A(i^ Hi) ^ t T ± 



xy{lh,l,k;mh,m-mh,mi) ^ (-i)' 



I' with A(l L' Ij) 
l+L'+Li even 



X y{l,lj,l';m,i - mh,mj,m, + rUj - mu) 

(2) 

(86) 



with 



7^, 



00 



X Rn, ij{r)r^ dr. 



(87) 

We can decompose 'R-n^i^^r.,u,n,i,,i,AkA) <M into a 
product of two independent one-dimensional integrals by 
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assuming the product ansatz 

gi{r,r') = gi^i{r) g2,i{r') (88) 

for the radial dependence (|83p . This simphfies our task 
to evaluate Eq. ([57]) greatly. It becomes 

'^nhlh.,nili,nj ^ '^iij'lj "^Uh lh.,ni li,l ' 

(89) 

with 

oo 

<'\,i,i'(fcA) = / J,{kAr)giAr)Rn,i,ir)r'dr (90) 



and 

oo 

<'Ln,ia = f Rn.ar')g2Ar')Rn^k(r')r"dr' . 



(91) 

VI. COMPUTATIONAL DETAILS 

All computations were carried out with MATHEMAT- 
ICA [g^I- In our essential-states model of Auger de- 
cay, the set of hole orbitals H comprises all five 3d or- 
bitals of krypton. The set of final states T consists 
of pairs of orbitals, the first is the 4s orbital and the 
second is a 4p orbital with magnetic quantum num- 
ber m 6 {-1,0,1}. The atomic orbital energies of kryp- 
ton are taken from Ref. [l^j: £34 = -70eV, £43 = -15eV, 
and e4p = -15eV. The effective charges (fM)) assume the 
values Z^d = 6.8, = 4.2, and = 4.2. The orbital en- 
ergies lead us to state energies (disregarding Hch + H^c) 
via Eqs. and This yields, for xuv pho- 

tons with LUx = 90 eV, a nominal photoelectron energy 
of np = Eo- +ujx = £3d + t^x = 20eV (H, [l^. We ob- 
tain a nominal Auger electron energy of SIa = E^ - E** = 
£As+£4p-£3d = 40 eV. The Auger decay width of a 3d hole 
in krypton is artificially set to Fbroad = 1-3 eV, which 
corresponds to a decay time of 500 as in accord with the 
data in Fig. 3 in Ref. [26]. We use this much shorter 
decay time to show the coherence in the laser-dressed 
Auger spectrum. The experimental value for the decay 
width is Fexpt = 88meV which corresponds to a decay 
time of 7.5 fs [l|]. 

With an approximation which we discuss below [see 
Eq. (|94p ] for the radial dependence (|88p. we determine 
the strengths of the dipole (|42p and Auger decay (|43p 
matrix elements. Using the decay width F|q^,=i from 
Eq. (|53)). we find from Eq. ([54]) the strengths Qv,Dxpt = 
1.10 with Foxpt and broad = 4.25 with Fbroad- The 
corresponding energy shifts follow from Eq. ([5^ . They 
are An^oxpt = -0.90 eV for (9„,cxpt and AR^broad = 
-13.53eV for Q^, broad- For a good agreement with 
the reference data [see Fig. [5] below], we employ the 
value Qv = 3.1 and the shift Ar = -0.68 eV. However, 
we set Ar = in all our computations because we have 
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FIG. 3: (Color online) Spherically integrated rms dipole ma- 
trix element d{kp) [Eq. (|93)) ] for the description of 3d pho- 
toionization of a krypton atom in the essential-states model 
of Sec. HVl 



chosen the orbital energies such that they correctly re- 
produce the (experimental) Auger and photoelectron en- 
ergies. The photoionization cross section of the krypton 
3d subshell for 20 eV photoelectron energy in Hartrcc- 
Fock-Slater approximation is read off of graph I in 
Ref. 'gg'I; it is about cthfs = l-SMbarn. The cross section 
without laser dressing a{uJx)\Q^=i is determined from 
Eq. dMl) using Fcxpt and letting Ar = 0. With Eq. ([70l). 
we obtain a dipole strength of Qd = 0.26. 

We use a xuv light pulse which has a Gaussian enve- 
lope with peak intensity Ix.o at i = and a full width at 
half maximum duration of tx = 500 as 

/x(i)=/x,oc-''"'(^^' . (92) 

The xuv electric-field envelope in Eq. (jlip follows 
from ex{t) = x/Stt a Ix{t)', its Fourier transform is also 
a Gaussian [56i] . 

VII. RESULTS AND DISCUSSION 

We devote this section to a computational study of 
our essential-states model and its parameters applied to 
M4^5iViiV2^3 Auger decay in krypton [1^53, [H, [3- It is 
motivated by a previous experiment [1J| which focused 
on the line group around 40 eV and related theoretical 
studies [25', [2^, [131 . We assess the accuracy of the approx- 
imations made in our essential-states model and compare 
with existing literature results. In future work (40| . we 
will explore laser-dressed Auger decay extensively. Ad- 
ditionally, we present the laser-dressed Auger spectrum 
for a much higher dressing-laser intensity than what has 
been used so far. 

To begin with, let us discuss the results for the rms ma- 
trix elements (|42l) and of the essential-states model. 
First, xuv absorption is determined by the rms dipole 
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matrix element in spherical polar coordinates (|42|) which 
is determined from the dh{kp,'dkp,'Pkr) [Eq. (|77p ] for 
all h € H. We take its modulus squared and integrate 
over the full solid angle to obtain the spherically inte- 
grated rms dipole matrix element, 



d{kp) 



J \d{h 

47r 



fcp y 



an 



(93) 



It is plotted in Fig. [3] After a steep rise at the edge (zero 
photoelectron momentum), it decays smoothly. Around 
the nominal photoelectron energy of Hp = 20 eV, the de- 
pendence of d{kp) on fcp is weak. Additionally, quantities 
such as the line shape function ([5^ decrease rapidly as 

soon as deviates appreciably from flp. 

Second, Auger decay is mediated by the rms two- 
electron matrix element (j43p . To construct it, we need 
the direct two-electron matrix element xj^{kx,'dkAiVkA) 
[Eq. ([86]) ] for all h eH. The rms matrix element in spher- 
ical polar coordinates is then denoted by v{kx, '&kA i "^^a )■ 
We examine three cases for the radial dependence (|83p 
and (|88p . This provides us with a good way to as- 
sess the quality of our approximation. First, following 
Refs. [2i,|23|, we use 



Second, we examine the reverse case, 



Qi{r, r') 



,,11+1 



(94) 



(95) 



and, finally, the exact case ([83]) . For comparison with 
Refs. [1^, m, [131 and for computational efficiency, we 
will use the crude approximation (I94p throughout. 

We display the rms two-electron matrix ele- 
ment C(fcA,0,0) in Fig. |4] with the viewing direction 
along the z axis which, in turn, is the direction of the 
linear xuv and laser polarization vectors. Our choice of 
direction agrees with Ref. [1^. The dependence around 
the nominal Auger energy of SIa = 40 eV — i.e., over the 
range plotted in Fig. [7] — is weak. 

To asses the impact of our omission of the two-electron 
exchange matrix element v,,z in the essential-states 
model of Sec. HVl we compute its rms value in spheri- 
cal polar coordinates vxikA^'&kAi VkA)- It is displayed in 
Fig. [5] along the z axis, vx{kA,0,0), for the approxima- 
tions in Eqs. (HH) and ((MI- The values of ux(fcA,0,0) 
in our momentum range of interest are roughly one or- 
der of magnitude smaller than the corresponding values 
of the direct matrix element at the same momentum in 
Fig. [1] for the approximation Eq. . The other case 
[Eq. ([55])] yields very large value for wx(fcA,0,0). Partic- 
ularly, these values are much larger than corresponding 
values for v{kA,0,0) which is unphysical. This compar- 
ison underscores that Eq. (|94p represents a reasonable 
approximation to the full Eq. (|83p while Eq. ([M]) does 
not. 
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FIG. 4: (Color online) The rms two-electron direct matrix 
element for krypton Sd-hole Auger decay ^(fcAjO, 0) of the 
essential-states model of Sec. IIVI Using the first radial de- 
pendence (|94p to find i;(fcA,0,0) yields the solid black line; 
the second radial dependence (|95|l gives the dashed red line; 
and the exact result is represented by the dotted green 
line. 
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FIG. 5: (Color online) The rms two-electron exchange matrix 
element i)x(fcA,0,0) for krypton 3d-hole Auger decay. See 
Fig- 111 for details. 
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2cJL 


0.00019 


0.0019 


0.019 
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TABLE I: Parameters for an 800 nm dressing laser. The 
vector potential amplitude Ai^ pip , the electric- field ampli- 
tude i?L,o P2p . the ponderomotive potential C/p (|37p. the 
magnitude of the maximum excursion ol psp , and the second 
argument in generalized Bessel functions for a broad range of 
intensities Jl o- 
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FIG. 6: (Color online) Auger electron spectrum (|63|l of 
laser-dressed krypton 3d-hole decay with an artificial width 
of Fbroad = 1.3 cV for a dressing-laser intensity of 5 x 10^^ 
The solid black line was determined using our theory; l5he 
dashed red line is the solid curve of Fig. 3 in Ref. [2^ scaled 
by a factor qIE\q = 2.9 x 10"* a.u. and shifted by -0.68 eV. 
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FIG. 7: (Color online) Auger electron spectrum ((63} of 
laser-dressed krypton 3d-hole decay for an artificial width 
of Fbroad = 1.3 cV and a dressing-laser intensity of 10^^ 
The solid black line was determined using generalized Bessel 
functions (I40() : the dashed red line was obtained neglecting [/p 
and using ordinary Bessel functions. 



We need to specify and characterize the xuv and op- 
tical light fields next. A present-day attosecond-pulse 
Ught source typically has a xuv peak intensity of at 
most /x,o = 10^^ at a photon energy of = 90 eV. 
The resulting ponderomotive potential ([57]) is 1.8 /zeV 
with a magnitude of maximum excursion (j38p of 8.2 fm. 
Clearly, the impact of the xuv field on photo- and Auger 
electrons can be omitted in excellent approximation. 
However, this approximation does not hold for the dress- 
ing laser — typically a Ti:sapphire laser system with near- 
infrared (nir) light of a wavelength of 800 nm and a pho- 
ton energy of wx = 1.55 eV — which delivers a large range 
of laser intensities. Exemplary data are given in Table [H 
For lower intensities — up to 10^^ — the ponderomo- 
tive potential (|37p and the vector potential amplitude 
[Eq. ([?!]) ] are small in relation to the laser photon energy. 
Also is negligible compared with the momentum of 
the photo- and Auger electrons. Therefore, we may ne- 
glect the influence of the ponderomotive potential [/p. 
Specifically, this amounts to replacing generalized by or- 
dinary Bessel functions in equations such as ((36)) . ((63|) . 
((S5|) . and ((75)) . For the Auger spectrum ((^ this was 
done in Refs. [2^, [l^] and was done also in this paper 
for laser intensities up to 10^^ This approximation 
is assessed by comparing Auger spectra from the expres- 
sion using generalized Bessel functions ()40p with spectra 
from the corresponding expression using ordinary Bessel 
functions. Excellent agreement is found. 

Finally, we are in the position to put together all in- 
gredients to compute the laser-dressed Auger electron 
spectrum ((63p in spherical polar coordinates of kryp- 
ton 3c? hole decay ^] V t\{kx,'&k^,^kj^^ integrated over 
the azimuth angle PaC^a) = Stt 7^A(fcA,0,0). As a veri- 
fication of our solution, we compute the Auger spectrum 
in Fig. [5] for a dressing-laser intensity of 5 x 10^^ as 



in Fig. 3 in Ref. )26ll . Our result agrees very well with 
the one in Ref. [26[ apart from an overall scaling fac- 
tor of E\ Q = 2.9 X ID"'* a.u. due to the xuv electric-field 
strength which was set to unity in Ref. [1^ . Finally, the 
spectrum in Fig. 3 of Ref. [2(| was determined with a 
nonzero value for Ar. We find that we need to shift it 
by a value of -0.68 eV [Sec. IVI) to achieve agreement. 
The structure in the figure is mostly due to the time 
dependence as the momentum dependence of the dipole 
and Auger decay matrix elements is weak. The differ- 
ences between both curves are ascribed to a somewhat 
different treatment of the matrix elements. 

We present the Auger electron spectrum for a dressing- 
laser intensity of lO" ^ in Fig. [7] [6i|. The spec- 
trum is first computed using generalized Bessel functions 
and a nonvanishing ponderomotive potential in Eq. ((63|) . 
Then, it is determined using ordinary Bessel functions 
and C/p = 0. There are clear differences between the 
spectra. The convergence with respect to the number 
of terms in Eq. ()40p is rapid; a summation from -2 to +2 
in Eq. ((10)) was sufficient We conclude that dressing-laser 
intensities around and above 10^'^ require an accurate 
treatment of the Volkov phase l(4l)) . We need 15 laser 
photon indices (Bessel functions) to account for absorp- 
tion and emission of laser photons in the sum ((61[) for the 
Auger electron spectrum. This is in good agreement with 
previous studies of iC-shell ionization of laser-dressed 
neon [2^ . argon [s^. and krypton [38!| atoms — however, 
with a very different theoretical approach — where 20, 12, 
and 5 photon blocks were required, respectively, to con- 
verge the calculations for the same dressing laser param- 
eters that are used here. Seemingly crucial for the neces- 
sary number of photon blocks are the decay widths of the 
inner-shell hole which were 0.27 eV, 0.66 eV, and 2.7 eV, 
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respectively. Our artificial value for the krypton 3d de- 
cay width of 1.3 eV lies between the decay widths of the 
krypton and the argon K shell vacancies. 

VIII. CONCLUSION 

We have devised and applied an ab initio theory for 
inner-shell xuv photoionization and subsequent Auger 
decay of laser-dressed atoms, a so-called two-color prob- 
lem. Our work aims at the study and control of elec- 
tron correlations — here manifested in terms of electronic 
decay — which is the most profound goal of attosecond 
science. The photo- and Auger electrons experienced 
an optical dressing laser which was considered to be 
intense but not strong enough to excite or ionize elec- 
trons in the atomic ground state. We used the Hartree- 
Fock-Slater (HFS) approximation as a starting point 
for the description of the atomic electronic structure. 
The HFS orbitals were then used to represent the full 
Hamiltonian. We employed a single configuration-state 
function to represent the ground state and singly- and 
doubly-excited states. The light fields were treated semi- 
classically and we use the strong-field approximation to 
treat the influence of the optical laser on the photo- 
and Auger electrons. The influence of the laser on the 
atomic ground-state electrons was neglected. The quan- 
tum dynamics of the problem was described in terms of 
equations of motion (EOMs). The EOMs were solved 
analytically for an essential-states model and a closed- 
form expression for the Auger electron amplitude was 
obtained. Furthermore, the xuv-absorption cross sec- 
tion of laser-dressed atoms and an expression for the 
laser-dressed photoelectron spectrum were derived. We 
applied our formalism to study the photoionization of 



a 3d orbital (M shell) of a krypton atom and its sub- 
sequent Mi^^NiN2^z Auger decay where the vacancy is 
filled with a 4s valence electron expelling a 4p valence 
electron. Following Ref . [2^ , we assumed an artificial de- 
cay width of krypton 3c? vacancies of 1.3 eV. The atomic 
orbitals were approximated by suitably scaled hydrogen 
wave functions circumventing the need for a HFS compu- 
tation. We discussed the approximations made and stud- 
ied the convergence of the Auger decay matrix element 
and the expansion in terms of generalized Bessel func- 
tions. We compared our laser-dressed Auger spectrum 
to literature results of Smirnova et al. [26| and found 
good agreement. Finally, we presented the Auger elec- 
tron spectrum for 10^^ NIR laser intensity. 

Our work opens up a multitude of future research per- 
spectives. We have devised a general ah initio framework 
which allows us to create simplified models of varying so- 
phistication tailored to model many physical situations in 
laser-dressed Auger decay. In this paper, we reduced our 
EOMs to an essential-states model which can be solved 
analytically and comprises sufficient details for a number 
of physical problems. In a forthcoming paper [40| . we will 
use it to investigate coherence and interference of Auger 
electrons and their control by a laser. In a next step, our 
model can be generalized to a few states, e.g., all states in 
a subshell can be considered, avoiding magnetic quantum 
number-averaged dipole and Auger matrix elements. 
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